Limit Theorems for Motions in a Flow with a Nonzero 

Drift. 



Albert Fannjiang 
Department of Mathematics, University of California, Davis 

Tomasz Komorowski 
Institute of Mathemathics, Maria Curie- Sklodowska University, Lublin 

February 1, 2008 

Abstract We establish diffusion and fractional Brownian motion approximations for 
motions in a Markovian Gaussian random field with a nonzero mean. 
Keywords turbulent diffusion, mixing, corrector. 

AMS subject classification Primary 60F05, 76F05, 76R50; Secondary 58F25. 
Abbreviated title Limit theorems for motions in a random field. 



1 



1 Introduction. 



The simplest model of a motion of a tracer particle in a turbulent medium, e.g. fluid, is 
given by the stochastic differential equation 



dx{t) = {v + eV{t, x)) dt + V2Kc/B(t), (1) 

where the particle is located at when t = 0. v & is a constant vector representing the 
mean velocity of the fluid, V(t, x) = (Viit, x), ■ ■ ■ , Vd(t, x)) describes the velocity fluctuations 
and is assumed to be a zero mean stationary random field. 

B(t) is a standard Brownian motion, whose presence accounts for the existence of the 
molecular diffusivity in the medium given by k > 0. The parameter e > is supposed to 
indicate smallness of the fluctuation amplitude in comparison with the mean drift. Turbulent 
flows of that nature have been extensively studied in Physics of Fluids starting with the work 
of the British physicist G. Taylor in the 1920's (see [0 and also P for more references on 
the subject). 

The question is to determine the long time behavior of the displacement x{t) from the 
statistics of the velocity and the size of the molecular diffusivity k. The effect of molecular 
diffusivity is in many cases negligible if k is small. We therefore set k = for simplicity of 
presentation. We assume also that the mean velocity v ^ 0. The case of vanishing mean 
velocity is presented by the authors in a separate paper (see @, 

More specifically we set 

yeit):=x(^)-v^. (2) 



We then have 

A theorem describing the limiting behavior of the particle, proven in various forms in 
||T0| , [12|, says that ^/^(t), t > 0, considered as continuous trajectory stochastic processes 
tend weakly, as e | 0, to a Brownian Motion whose covariance matrix is given by the so-called 
Taylor - Kuho formula 

oo 

D*^ = J {E[Vi{t,vt)V,{0,O)]+Em,vt)Vj{0,O)]}dt (3) 



(cf. 0). 

The version of the above theorem proven in [IC] for v = 0, requires, among others, the 
strong mixing assumption for the velocity field in the temporal variable. The key point of 
the argument is an approximation of the forward particle location valid for sufficiently large 
times that only takes into account the information available prior to a given instant, the 
so-called frozen path approximation. The temporal mixing property of the field guarantees 
then that the future particle displacement is almost independent of the past information 
thus, in consequence the diffusion approximation holds. This type of theorem can be also 
shown for a field with v ^ provided that its fiuctuation is strongly mixing in space (see 0). 
In this case one can project the spatial domain into time domain via the change of variables 
X := X — vt, cf. (H), and then, using the fact that the mean drift size dominates, for small 
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e, the amplitude of fluctuations it is possible to translate the spatial mixing property into 
temporal one for the velocity field considered in the new coordinates and, as a result, the 



argument of |]10[ can be applied. 

In this paper we set out to give a description of the limiting behavior of scaled trajectories 
for a class of Gaussian fluctuation fields with power-law spectrum. This family of fields 
appears frequently in the mathematical theory of turbulence and is motivated by the seminal 



work of Kolmogorov |TT|, see also 0] for the shear layer flow case. In many interesting cases 
the fluctuation fields with that kind of spectrum do not satisfy the mixing condition either 
in time or space, thus the existing versions of limit theorems are not applicable. 

In what follows we assume that V(t,a;), {t,x) E R x W^, is a d-dimensional, time- 
stationary, space-homogeneous, time-Markovian, isotropic Gaussian velocity field defined 
over a probability space {Q, V, P) whose two-point covariance matrix R = [Rij] is given by 
the Fourier transform: 

Ri^{t,x)=E[V,{t,x)Vj{0,0)] = J Co(k-a3)e-l''l"'*A,,(k)rfk. (4) 

n't 

Here and in the sequel we let Co(0) := cos0, Ci(0) := sin0. E stands for the mathematical 
expectation calculated with respect to P. The spatial spectral density of the field R = [Rij] 
is given by 

where a : [0, -|-oo) — >• the so-called ultraviolet cut-off, is a certain compactly supported 
continuous function. The factor 6ij — fcjfcj/|kp in (|^) ensures that the velocity field is 
incompressible, i.e. V- V(t, x) := diVi{t, x) = 0. The function exp(— |kp^t) in @ is called 
the time correlation function of the velocity V. The spectral density R(k) is integrable over 
k for a < 1. The ultraviolet cut-off is then needed to avoid divergence of the integral over 
|k|. The parameter a is directly related to the decay exponent of R. Namely, R(a;) ~ |a:;|"~^ 
for I a; I ^ 1. As a increases to one, the spatial decay exponent of R decreases to zero and in 
consequence the strength of spatial correlation of modes increases. For a > 1, the spectral 
density ceases to be integrable and needs an infrared cut-off, i.e. the origin should lie outside 
of the support of a. On the other hand, in order to apply the results of or [|l^ one needs 



to assume that (3 < 0. Otherwise, i.e. when (3 > 0, the velocity field lacks the spectral 
gap and thus the strong mixing property (cf. In light of the foregoing discussion we 

therefore restrict our attention to the case a < 1, /? > 0, which corresponds to the velocity 
field having arbitrarily long scales but not the strong mixing property. 

In addition, the particular form of the spectrum of the field allows us to assume without 
any loss of generality that it is jointly continuous in both t and x and is of C°° class in x 
a.s. (see e.g. 0). 

One can then pose the following question: what is the region in the [a, (3) plane where 
the classical turbulent diffusion theorem, with the Taylor - Kubo formula holds? It is 
easy to find the necessary condition by imposing the convergence of the integral appearing 
in d). We have 



D*^ = Jp^.it, vt) dt = J U, - J expHM'H)co{kvt) dt rfk, (6) 



Rd . . / . . 
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where := k- v. After a straightforward calculation of the integral over t one deduces that 
the rightmost part of (|) equals 

.(Ik|)|k|- (^_^ _ ^ 



|k|2a+d-2Qk|4/3 + 4) 



Elementary calculations (cf. the computations made after and (^Bj) below) show that 
the finiteness of the expression (0) leads to the conditions 

a + p <1 when p < 1/2 (8) 



or 

a<l/2< p <1. (9) 

It turns out that these conditions are also sufficient to claim that the particle trajectory 
approximates the Brownian motion in the sense of weak convergence of continuous path 
stochastic processes. This is essentially the contents of our Theorem |I] below. 

One can also ask what happens in the case when conditions {^) are violated. The fact 
that the integral in (P) diverges suggests "slower" than scaling of the temporal variable 
in (0). We define 

yM--=^{^)-v^ (10) 

for a certain suitably chosen 5 G (1/2, 1). The limit of the stochastic processes given by the 
trajectories turns out then to be a fractional Brownian motion. This result is contained in 
Theorem |^. 

It is interesting, in our judgement, to compare the results of the present article with those 
of and 0. It was proven there that the diffusive regime for the motion in the field without 
the drift, i.e. v = 0, comprises of (a,/3) satisfying a + P < 1 (see 0). When, on the other 
hand, a + P > 1, the properly scaled trajectories approximate a fractional Brownian motion 
(see 0). As a result we see that inthe present setting the diffusivity regime is enlarged 
for P < 1/2. We call this phenomenon, after [Q, the sweeping effect. It can be explained 
intuitively as follows. For a given wavenumber of magnitude |k| one can see from (D that 
the relaxation time needed for a significant temporal decorrelation of the fiuctuation field 
is of order r ~ |k|^^^. On the other hand the sweeping time necessary for a significant 
decorrelation related to the mean drift convection is of order of magnitude tv ~ |k|~^- It is 
much shorter than the relaxation time when /3 > 1/2. In that regime therefore the mean drift 
is responsible for the chaotic behavior of the velocity field and in consequence the diffusive 
approximation holds when a < 1/2. On the other hand for P < 1/2 the effect of the mean 
drift is negligible and the limiting particle motion is the result of a subtle balance between 
the spatial and temporal chaotic properties of the field as described in 0, |§] . 

To simplify the presentation we focus our attention only on the case when a, P satisfy in 
addition 

a + 2p<2. (11) 

This condition can be proven inessential for the validity of our results. The more technically 
involved case when a + 2/3 > 2 will be the subject of a subsequent article. 
We summarize our results in the form of the following two theorems. 
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Theorem 1 Suppose that 

1) vj^O and (|n|) holds. 

2) Y{t,x), {t,x) G i? X R^, is a zero mean, Gaussian velocity field whose correlation is 
given by (|) and (||). 

Then, for (a,/?) as specified by conditions (||) and (|) the scaled trajectories given by ^ 
converge weakly, when e I 0, as continuous trajectory stochastic processes to a Brownian 
motion whose covariance matrix Dq is given by @. 



Theorem 2 Assume conditions 1) and 2) of Theorem^. Then for 1 > a, [3 > such that 

a + l3>l when 1/2 < a <1 (12) 



one can find a unique Sa^p G (1/2, 1) and H 



25. 



'a. 13 



such that the scaled trajectories given 



by ([ToD converge weakly to a fractional Brownian Motion B/^(t), t > 0, i.e. the unique 
H -self-similar, Gaussian process with stationary increments (see ^\). We have 



with 



where r^,/? equals 



D 



E[B^(t)®B^(t)] = D„,^t2H^ 

k® k\ a(0)rfk 
IkP J |k|^-i ■ 



afi 



ra,/3 I 



(13) 
(14) 





- 1 + 


k 


2/3 




]^ 2a+4/3-l 



when /? < 1/2, 



°^ '"^ when (3 > -, 



1.2 \U\2a-l 

(|k|2 - kl) (e-^'^koikv) - 1 + |k|) - 2|k|A;^ (e-^'^kiikv) - kv 



In addition, d^n equals 



|k|2 + 4)|k|2«-l 



when 3 = -. 
^ 2 



a + 2(3-1 
1 

2a 



when (3 <l/2, 
when 13 > 1/2. 



(15) 
(16) 



2 Preliminaries. 



The Spectral Theorem for real vector valued random fields (cf. e.g. |]l|) implies that there 
exist two independent, identically distributed, real vector valued Gaussian spectral measures 
V^t, ■) = (V^,i(t, ■),■■■, Vi^d{t, ■)), ^ = 0, 1, such that 



\{t,x) = / Voit,x,dk), 



(17) 
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where 

Vo(t, a;, dk) := co(k ■ x)Vo{t, dk) + ci(k • x)Vi{t, dk). 

In the sequel (cf. e.g. Lemma |lD we also deal with the spectral measure given by 

Vi(t, X, dk) := -ci(k • a;)Vo(t, dk) + co(k ■ £c)Vi(t, dk). 

The field V is Markovian in the following sense. For any function ip G S{R'^,R) and 
i E {1, ■ ■ ■ ,d} we have 



E 



I ^P{k)Vi,it, X, dk) I V_oo,J = / e-l'^l"'(*-^)^(k)V^,(s, a;, dk), 1 = 0, 1. (18) 



Here, for any a < b, Va,b denotes the cr-algebra generated by V(t, x) with (t, as) G [a, 6] x i?'^ 
and Vi^i{s,x,dk) is the relevant component of Vi(s, a;, rfk). In what follows we also write 
L2,for L'{n,Va,b,P). 

We now introduce the key concepts appearing in the proofs of Theorems |I| and namely 
Xa, the X-corrector and Ua, the X-convector. We set 



(t, y) := J e-^(^-*)E [V(s, -J^s + ?/) | V-oo,t] ds, (19) 



UA(t, t/) := V(t, + y) ■ VxA(t, y). (20) 

The differentiation used in ([20| ) is understood in the mean square sense. 

The procedure defining x\ and Ua requires some justification at least to guarantee 
smoothness of the corrector field appearing on the right hand side of (pO]). This techni- 
cal point is explained by the following. 

Lemma 1 For arbitrary A > the fields xx, Ua defined above are integrable for any 
p > 1. LP norms of those fields have the following hypercontractivity property: 

iE\xxit,y)\n'^'<C(B\xxit,y)\'y^\ (21) 



(E|U,(t,y)n^/^ <C{B\V,{t,y)\'y^\ (22) 

where the constant C > does not depend on X. 

One can also select modifications of x\ o-nd \J\ which are jointly stationary in the strict 
sense, continuous in both {t,x) and C°° in x. 

The spectral representation of xx is given by 

Xx{t,y)=xf\t,y) + x['\t,y) 

with 

xf{t.y)-= lQik,X)V,{t,vt + y,dk)t = l,2 (23) 

when 

and ^ 

C'2(k,A) := (|k|2/3 + A)2 + 4 ^2^) 



6 



The proof of the lemma is fairly standard. We sketch here only its main points referring 
the reader interested in details to the relevant literature. integrability and inequalities 
(PI), (l2^ ) for the fields in question follow from the hypercontractivity property of norms 
in Gaussian measure spaces related to the velocity field V (see e.g. Theorem 5.1 and its 
corollaries in p). 

The existence of regular versions of the fields follows from the fact that there exists h > 
such that for any integer > we can find a constant C^^h > depending on h, N only for 
which 

E|D-V(t, x) - D^Y{s, y)\^ < Cn^ [\t - s\'' + \x - y\' 

\in\=N 

for all [t, x), (s, y) & R X R^. Here for any integral multiindex m = (mi, ■ ■ ■ , m^) we define 
jjm ._ gmi ...gmd_ ^j^Q pg^j^ gj^^ ^]^g modifications of the fields with required degree of 
regularity using the aforementioned Theorem 5.1 of ibid in conjunction with Kolmogorov's 
classical result on the existence of continuous trajectory modification of a stochastic process 
(the random field version of that criterion can be found in e.g. as the Corollary to 
Theorem 3.2.5). 

Joint stationarity of the fields in question can easily be verified by an application of the 



results of [15 



3 Proof of Theorem 

We define a scaled corrector along path by 



Xe{t):=exe^{^,ye{t)) (26) 



and a scaled convector along path by 



U,(t):=U,.(3,i/,(t)). (27) 



Thanks to the divergence free structure of the velocity field V the processes Xeit), t > 0, 
and Ue(t), t > 0, are jointly strictly stationary (see Theorem 2, p. 500, of |T^). 

The importance of the concepts of a corrector and convector is highlighted by the follow- 
ing lemma. 

Lemma 2 We have 

t t 

Veit) = -Xe{t) + Xe{0) + J Xe{s) ds + J U,(s) ds + M,{t). (28) 



Here Mi,{t), t > 0, is a Brownian motion whose covariance matrix is given by where 



In addition, there exists 7 > such that 



lim£-^E|xe(t)r = and lim£-^E|Ue(t) T = 0. (30) 
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Proof. First we show (28). We start with the calculation of the so-called pseudogenerator 
of the path corrector, i.e. the process given by (cf. Chapter 3 of p 



^ ^ 510 



E 
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(31) 



where the limit is understood in the sense. 

It is easy to observe via Taylor's expansion of Xei^i v) about y^{t) that 



Xe{t + 5) = ex, 



't + 5 



£ 1 9 



,y,it)\+Rs{t), 



(32) 



where 



Here o{6) is understood as a function with takes values in the space random vectors with L'^ 
integrable components for an arbitrarily chosen p > 1 and such that o{6) / 6 vanishes in the 
L'P sense as 5 | 0. 

Conditioning the first term on the right hand side with respect to the cr-algebra V_oo,t/£2 
we deduce from (|T9|) and (Eol) that 



eE 



ft + 6 

Xe^ ^-,2/e(i) I V_oo,t/£2 



(33) 

The last expression, after applying Taylor's expansion this time about t/e"^, can be rewritten 

as 



Xe{t)+S 



Xe{t)--V[-,,v- + y^it) 



6 ve^' e"^ 



+ 0(5). 



In consequence, 



Lxeit) = Xe{t) - -V v\ + y^{t) ) + U,(t). 



Using the results of Chapter 3 of [T^ we obtain that 

t 

M,{t) = Xeit) - j Lxeis) ds 



(34) 



(35) 



(36) 



is a continuous trajectory vector valued martingale. We write = (Mg^i, ■ ■ • , M^^d); a 
similar convention will be used for other vector valued fields appearing with subscript e. 

We notice that thanks to the aforementioned results of |T^ the joint quadratic variation 
of martingales M^^i(t), M^j{t), t >0, i,j = 1, - ■ ■ ,d, equals 



(M,„M,,,)i = J L{M,,{s)M,^,{s)}ds. 



(37) 
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After elementary calculations we get 



L{M,,,(t)M,j(t)} = lim 
E ' 



E 



lim 

<5iO 



HO 6 



(3J 



5 



-[x.,.(t)i:x.j(t) + x.,,(t)i:xs,^(t)]. 



Using expansion (|3^ to represent ^^(t + 5) in ( pSj) we obtain that the utmost right hand 
side of (|38|) equals 



L.^jit- y,{t)) + f/,,i(t)x.,,(t) + XeAt)Ue.,{t) - XeAi)LXeAt) " , (39) 



where 



E 



Xe,ii. ^2 ; 



Xe,i{^,y)Xe,j{^,y) 



The corrector fields appearing in ( ^OD are zero mean Gaussian therefore 



(40) 



E 



(ill \ 



Em„- + E[x.,.(0)x.,,(0)] 



where 



= e^E 



E 



X. 



Xe,i(t)+Lj(t) - -V^, f4,'"4 + Veit) 



(41) 



With the help of (^) and (Q) we deduce that 

Le,ij{t; y,{t)) = Be,i,,{t) - E5,,,,,(t) 

where 

Be,i,j{t) ■-- 

Combining (^I]) with (|39| ) we deduce, using (|35|) , that 

L {M,,i(t)M,,,(t)} = -E5,,,j(t) = -EBe.jiO), (42) 

with the last equality following from stationarity of the relevant processes along the particle 
trajectory that holds due to the results of [jl5|. A direct computation, using the formula for 
the expectation of a product of Gaussian random variables, shows that the right hand side 
of ( ^2D equals D^^ij, where := [-De,^^] is given by (p9D. This concludes the proof of (pSj). 

Now we prove formula (|30|) . Without any loss of generality we assume that v = 
(l,0,---,0). 

We then have 



B\xe{t)f = E\xeiO)f<Ce' I 



dk 



\k\<K 



[(£2 + |k|2/3)2 + A;2] |k|2"+rf-2- 



(43) 
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If we represent k = {ki, 1), where 1 G i?'^ ^ the rightmost part of (43) can be rewritten as 

dkidl 



^kl+\\\^<K 

<C' 



{[e^ + {kl + |l|2)/5]2 + kl} {kl + |l|2)"+(<i~2)/2 

e'^l'^-'^dkidl 



(44) 



{[£2 + [kl + /2)/3]2 + kl} {kl + PY+(d-2)/2 ' 

After changing variables in the last integral according to the rule / = r cos (p, ki = r sin (p we 
infer that the leftmost part of ( ^5[ ) can be estimated from above by 

C 



e^drdcp 







[(^^2 _j_ j,2/3-)2 _|_ ^202j^2a-l ' 



Here C denotes a generic constant independent of e. 

After an elementary calculation we find that the last integral equals 



K 



arctan 



7rr 



dr. 



For a, 13 as specified by conditions (| 
that this expression vanishes as e | 0. 
To show (|30|) we first observe that 



r2"(52 + r2/3) \^2(£2 +r2/3)^ 

I, (^ the Dominated Convergence Theorem implies 



E |Ua(0, 0)1' < CE |V(0, 0)1' E |Vxa(0, 0)|' 



The gradient of the corrector satisfies 
E|Vxa(0,0)|'< 





k 


2 


[(A + 


k|2/3)2 


+ kl] 


k 


2a+d-2 



dk. 



(45) 



(46) 



Estimating precisely in the same way as we did for (^31) we infer that the right hand side of 
(H3) is less than or equal to a constant times 



K 



r2a-2(^X _^ ^2/3) 



arctan 



vrr 



^2(A + r2/3)^ 



dr. 



This expression remains bounded as A | if only a + 2/? < 2. Consequently, 

lim sup E |Ua(0, 0) |' < +oo. (47) 

AiO 

(^) holds also for any p-th absolute moment of Ua with p > thanks to Lemma |l|. 
We also have 



E 
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t s 



(4J 







t s s 

+ -JdsJ ds, I E|(v(|,^| + y,(.2)) -VU,. ^^^,y^{s,)^yv,. 



Si 

The first term of (HHf) equals 



2jdsjE |U,2 





(49) 



due to tlie stationarity of the path convector. 

Thanks to the incompressibihty of the velocity field we have E\Jx{t, y) = 0. That and well 
known properties of conditional expectations of second degree polynomial-like functionals of 
a Gaussian field then imply that 



Vx{t, y- 0) = V(t, vt + y; 0) • Vxx{t, y; 0). 
Here V(-,-;0), Vxa(-,-;0) denote the orthogonal projections of V(-, ■) and Vxa( 

-'^-oo.O- 

In consequence, ( ^91) equals 

t s 

2 1 rfs I E I v'^; o) ■ Vxe^ 0; o)) ■ U,2(0, 0)} ds,. 



(50) 
onto 



(51) 







Using the Cauchy Schwarz inequality we estimate (|5TD by 



C,{E\VA0,O)\Y\E\V{0,O)\y/Ujds J E 



1/2 



^0 



Vx,2 (^,0;0 



dsi 



(52) 



The fourth moment of remains bounded as e | 0, see (^71). Therefore, the expression 
in (|52[) can be estimated by a constant times the last factor of ([5^). The latter, with the 
help of (|TB|) and (||), equals 



\k\KK 



e^[l-exp(-«^)] 

|k|2/3 



X 



|k|2a+d-d [(^2 + |k|2/3)2 + ^2] 



(ik (is 



1/2 



(53) 



The expression (^3|) tends to by virtue of the same change of variables as the one used 
above and the Dominated Convergence Theorem. 

Finally, we consider the second term of (0). After a simple change of variables it can be 
rewritten as 

t s s 

-JdsJ ds, |E|(v(^,i;^ + y,(si)) • VU,2 (^£^,^^(,0;^)) ■U,2(0,0)}ds2. (54) 

si 
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An analogous argument to the one we used to deal with (^91) shows that (54) vanishes as 
e I 0, provided that 



4e 



-si 



VU,2 i%n-o]ds 



(55) 



<4e 



S — Sl 



+ -E 



tends to zero as e | 0. Here 



VU<?(i|,OiO) :=V 



v(3.t.5;0).Vx-«(f|.OiO) 



i = 1,2. 



(56) 



Using spectral representations (p3D, (^), (^) and (0) we can write the p, q-th entry of the 
matrix defined by (|56D as being equal to 



exp 



S2 



-/expj- 



S2 



Ci{'k,e )kVi^i^p [0,v—,dk) ■ / /^exp 



S2 



a(k,£^)kA;,\/,,JO,t;^,rfk ■ / exp 



|l|2/3 



S2 



^2 I Vi (^0,1;^, rfl 



Vo(0,«^,(il 



Hence, each term of the rightmost part of (|55D equals 



exp 



xCi(k, £^)Ci(k', e^)W{dk, d\, dk', dl'; S2, 4) 
with W the signed measure given by 



d 

EE 

p,q=l 



S2 



S2 



i-iykgVi.,,p 0, v^, dk k ■ Vi 0, t;^, rfl - k^V,,p 0, -j;^, dk k • Vq 0, i;^, d\ 



S2 



S2 



iyk'^Vi-i,p I 0, v^, dk' 1 k' ■ Vi ( 0, v^, d\' 1 - I 0, v^, dk' I k' • Vo ( 0, t; J, d\' 



The above expectation is to be calculated treating the spectral measures as formal Gaus- 
sian random variables and subsequently applying the rules of calculating the expectation of 
products of Gaussians with the following relations for the covariance 



E 



V,(0, X, dk) O V,/(0, x', dk')] = 5i,rCo(k ■ {x - a;'))R(k)5(k - k')rfkdk'. 



After a straightforward computation we obtain that 



W{dk, d\, dk', dV; S2, 4) = 2\k\'^co{k_ 



S2 - S2 



)co(/: 



S2 - S2 



)trR(k)k'^R(l)k 
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x5{k - k')5{\ - r)dkd\dk'dV. 
Thus each term on the utmost left hand side of (^5]) can be bounded from above by 

C / / CfiKe^) , niL^nnnmU.,. dl (57) 



(|k|2/3+ |l|2/3)2(|k||l|)2a+<i-2 
|k|,|l|<K 

with the constant C independent of e. For a + 2/3 < 2 (0) tends to zero, as £ J. 0, by virtue 
of the Dominated Convergence Theorem and the proof of the Lemma ^ is finishedn 
The proof of Theorem^. By virtue of we write 

y^{t)=M,{t)+R,{t) (58) 

where the remainder term Rs{t) consists of a stationary, i.e. — Xe(^)) and an additive part 



RaAt) ■= J Xe{s)ds + J lJe{s)ds. 



According to Lemma ^ Ms{t), t > 0, is a Brownian motion whose covariance matrix, 
given by (p9|), converges to Dq as e | 0. In consequence, M^^t), t > 0, converges weakly, 
when e I 0, as processes with continuous trajectories (cf. e.g. 0). 

In order to conclude the proof of Theorem |l| it remains to be seen that for arbitrary 
T > 0, 5 > we have 

limP[ sup \Reit)\ >5]=0. 

'^iO 0<i<T 

The assertions of Lemma ^ imply that we only need to show 

limP[sup \Xe{t)\>S] = 0. 

£-10 0<t<T 

Since Xe is stationary we can write 

P[snp \xe{t)\>6]<^P[ sup \Xe{t)\>S]. (59) 

0<i<T E-" 0<t<Te2 

Using again (^) to represent Xe we can estimate the right hand side of (|59D by the sum 

sup [\R^^^(t)\ + \yM]>S] + ^P[ sup \MM>S] + ^P[\Xem>5]. (60) 

£ 0<t<Te2 £^ 0<t<Te^ £ 



The last term in ( pH ) can be estimated using Chebychev's inequality by 

nxe{o)n 



(61) 



Choosing p so that 'j -p > 4 with 7 given by Lemma |^, we see that the expression (^) tends 
to as e J. by virtue of (0). 
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The second term of (|60|) vanishes as e | thanks to 



lime-'E\M,{e'T)\' = 

ej,0 



and an elementary martingale inequality 

T 

0<t<T£2 

The first term of (^) can be estimated by 



;P[ sup |M,(t)|]>5]< 



C^<!^EMO)p + ^E|U,(0)p + |p 



E\Me{e^T)\^T 



1 

J \V{s,vs + y{s))\ds > 



(62) 



(63) 



The first two expressions in ( |63D tend to thanks to Lemma ^ The last term of (|63| ) vanishes 
as £ 1 0, provided that p > 2. This is due to the fact that 



lim sup E 



1 

J |V(s, vs + y{s))\ds 



< E|V(0,O)|^rP < +00. 



This finishes the proof of Theorem |l| □ 



4 Proof of Theorem @. 

Throughout this section we assume that the scaled trajectory is given by (|TU]) with 6a, /3 
given by (|T5|), ([T6|). To avoid cumbersome notation we suppress writing the subscript of 6. 
Define the path corrector and convector by 



and 



t 



U,(t) ■.= \],.s{—^,yM- 
The following lemma can be proven exactly the same way as Lemma |^. 

Lemma 3 The scaled trajectory of a particle satisfies 

t t 

y^it) = -xe{t) + Xe(0) + / Xe{s)ds + E^^^-'^ J lJ,{s)ds + £^-'M,(t), (64) 



where M^{t) is exactly as in the statement of Lemma 

Remark. Calculations analogous to those made after (^) show that for f3 < 1/2 the 
exponent 6 is chosen to make E\e^~^ M^it)]"^ of order 0(1). Indeed, 



E\e^-^M,{t)\'' = e''~'''ti-D,t ^ e 



2-284. 



^2(l+5i^) 



(65) 
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and 



In the case when /3 > 1/2 we can deduce by the same argument that the terms on the left 
hand sides of (pSf ) and ( |66D vanish as e J. 0. 

We note here that 6 < lforQ! + /3> 1 and j3 > 0. Thanks to the assumption that 
a + 2/? < 2 we have (cf. (|6|)) 



£^-^^E|Vx.2.(0, 0)p ~ E^-^' i as £ i 0. (67) 



This together with (^) imphes that the term on the right hand side of (|6^) involving the 
convector vanishes in probability as e J, 0. In addition, we deduce by ( |67D that 



X> 



t 



Here and in the sequel we use the generic notation Oeit) for any stochastic process which 
satisfies limE sup |o£(t)P = 0, for any T > 0. 

With this, (|6^ ) becomes 

i 

Veit) = -exe^^ ( J^, o) + eXe^^iO, 0)+eJ Xe^^ o) ds + e'"'M,{t) + o,(t). (69) 



Lemma 4 We have the following expansion of the corrector: 



where is a Brownian motion which satisfies 

i) the i,j-th entry of its covariance matrix is given by (|^). 

ii) there exists 7 > such that 

lim£2(i-5)-7E M,{t) - M,{t) ^ = 0. (71) 

ej.0 

Proof. The proofs of ( [TDD and part i) go along the lines of the proofs of the corresponding 
parts of Lemma so we leave this argument out. Actually, the fact that the spatial argument 
is fixed significantly simplifies the proof in the present case. 

We focus on the proof of part ii). We first calculate the pseudogenerator 

C{t) := £2(1-5)/: M,{t) - M,{t) ' 
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In light of (^) and (|64D it is clear that 

m = e'C\Zit, y^m' - e'Zit, yM " ^Z(t, yM (72) 

where ^ ^ 

Z(t, y) := x,2* 2/) - X.^^ o) . 

The first term on the right hand side of ( [7^ ) can be computed as limit of 

E{|Z(t + r,?/,(t))|2 I V_^,/,..} - |Z(t,t/,(t))p 

e lim — . 

rio r 

A direct computation similar to that performed after (|3T|) shows that 



m = 2e'^'-'^ / ' L/,^ , ,.,, [1 - co(k . xM)] dk. (73) 





k 


2/3a(|k|) 


lk 




25 _|_ 


k|2/3)2 + A;2] 



Here, recall, cc£(t) = vt + y^(t). For 2/3 > 1 computations analogous to those performed 
after (|i3| ) or ( ^6]) convince us that the right hand side of ([73| ) vanishes as e | 0. 

To estimate the mathematical expectation of the right hand side of (|73D when 2/3 < 1 we 
break it into two disjoint integrals: the first over |k| < e^^^f^, with < g < 1, and the second 
over its complement. 

The expectation of the second integral can be estimated by 

» |K.|2/3+2-2a-d 

2^'^'"'^ / w 2 II 12^2 ,21 dk. (74) 

K>\k\>e«e/I3 

Employing again the procedure used to estimate the last part in ( |i5| ) we infer that ( [7^ ) is 
less than or equal to a constant times 

^ ^2/3+1 

e2(l-<5) / !^ 

7 j,2a(rW _|_ £25-12 ' 

eT 

With the change of variables r' := e^^r we conclude that the above expression vanishes, as 
e I 0, at the rate 0{e^) for some 7 > since Q < q <1. 

The expectation of the first integral can be estimated as follows: 



t ^ 

- sin(k • Xe{s))k ■ V f-^, Xe{s] 





k 


''<\k\) 


|k 




25 + |k|2/3)2 + ^2] 



|k|<e^e//5 



dk 



<Ce^-^' / j ' ,^ — -2Tc/k. (75) 
J \ie^^ + k 2/3)2 + A;? 

Via the same type of argument as used to estimate the previous integral we find that 
]) is bounded from above by 

^ ^2/3+2 



' j,2a(^j,2f3 _|_ X)2 



16 



The divergence of the integral in (ffSD, if occurs, can be made arbitrarily slow by choosing q 
sufficiently close to one. Thus the expression in ( [75| ) vanishes for a + 2/5 < 2. This concludes 
the proof of Lemma ^ 

As a consequence of Lemma ^ we obtain 

x,{t) = y,{t) + Oe{t), (76) 

where 

yM=^ r V(^>0)^^ (77) 

JO 

is a Gaussian process whose covariance matrix equals 

E [y^{t) ® y^{t)] = 2£^ J J ^(^'' '"^') (78) 



An elementary calculation shows that the right hand side of ([78| ) tends to T)t as e | with 
D given by (p!^ . Hence y^ converges weakly as e | to a fractional Brownian motion Bnit) 
whose covariance is given by (|TB|). 
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